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Conference summary: Mass loss from stellar clusters 
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1. Introduction: topics on this conference 

This conference dealt with the mass loss from stars and from stellar clusters. 
In this summary of the cluster section of the conference, I higlilight some of 
the results on the formation and the fundamental properties of star clusters 
(§ E^)) the early stages of their evolution (§ [3]]) and go into more detail on the 
subsequent mass evolution of clusters (§[!]]). A discussion on how this may, or 
may not, depend on mass is given in § [5?] Obviously, there will be a bias towards 
the topics where Henny Lamers has contributed. Some of the contributions to 
these proceedings have already reviewed extensively the topics of clusters mass 
loss and disruption, so I will try to fit these in a general framework as much as 
possible. 



2. Cluster formation and fundamental properties 

Star clusters are considered to be the lowest step in the hierarchical star forma- 
tion scenario (see Elmegreen in these proceedings). This hierarchy is thought to 
result from a combination of (gravitational) fragmentation and turbulent com- 
pression in the interstellar medium (ISM), resultin g in scale-free pow er spectra. 
These have been observed for the Milky Way ISM (Dickey et al.l 200ll ) and seem 
also to exist for y oung star fields in M33 (lElmegreen et al" 2003 ) and NGC 628 
(jElmegreen et al.li2006. ). Fragmentation in a turbulent medium yields power- 

jutions. Powe r-law 
Salpetei"l955') and 



law distributions at all scales for the mass and radius distri 
mass functions are found for the stellar ma ss function (e.g. 



for the cluster initial mass function (CIM F) (jBattinelli et al. 



1994. Zhang fc Fall 



1999l: lde Griis et al.ll2003l:lBik et al.ll2003l1. with some evidence for an upper-mass 
truncation ( Gieles et al.l 2006a b; Powell et al. 2006 ^ and Haas et al. in these pro- 
ceeding s). A truncation is also found for the mass function of giant molecular 
clouds ( Williams fc McKe3ll997l : iRosolowskvl |2005| ) . but a sohd theoretical un- 
derstanding of the top end of t he CIMF is sti l l lack ing. When extrapolating the 
cluster mass function to = l. lde Wit et"all (|2005l ') show that 4% of the O stars 
should form in isolation, which is in agreement with their observations. On a 
larger scale, complexes of stars and clusters are als o characterized by scale-free 
power -law distributions for luminosity and radius (Ivanov 20051 : Bastian et al 
20053). 



The radius distri bution of clusters in M51 is well described by a power-law 
( Bastian et al. 2005bl ). For the interstellair gas clouds there is a clear relation 
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between radius and mass of the form M oc 



complexes (jElmegreen et alj 120011 : iBastian et al. 



which i s found back for stellar 
2005al ). This relation breaks 



down at the level of individual clusters, for which observatio ns are more con- 
sistent with no relation between mass and radius at all (e.g. IZepf et al]|l999l : 
iLarsen' '2 004). For the ern bedded clusters in the Milky Way this relation is still 
present (Allen et al.ll2007l ). suggesting that the relation gets wiped out by a com- 
bination of effects in the early stages of evolution (e.g. gas removal, dynamical 
mixing, tidal truncation etc.). However, the range of masses of these clusters 
span only a few orders of magnitude, which makes it dangerous to extrapolate 
this result to all masses. Several scenarios have been introduced to explain the 
scatter relation between cluster radius and ma ss, involving a mass dependent 
star formation efficienc y (Ashman &: Zepf 200lh o r the early evolution were gas 
expulsion plays a role ( Goodwin &: Bastian 20061 ). Still, a concluding explana- 
tion can not be given yet. 

Star clusters with masses in excess of a few times 10^ Mq tend to be slightly 
larger than the typical radii found for young clusters and globular clusters 
(~ 3pc). The mass-radius relation for these super massive clusters is well de- 
scribed by the Faber- Jackson relation for elliptical galaxies, when e xtrapolated 
down to these masses ( Hagegan et al.|[2b05l : iKissler-Patig et 311120061 and in these 
proceedings). When clusters form through merging of multiple clusters into a 
single massive object, the final rad ius is somewhat larger tha n the typical radii 
(~ 3pc) of clusters of lower mass ( Fellhauer Kroupa l2005l ). This points at a 
difference between the formatio n process of clusters with M < 10^ — lO"^ Mq and 
those with M > 10^ - 10^ Mq (iKissler-Patig et al.l[2006l l. 



3. Early evolution 



The early evolution (first ~ lOMyr) of clusters is dictated predominantly by 
the removal of left-over gas (not used for star formation) by stellar winds from 
early- type stars. This removes a significant fraction of the binding energy of the 
embedded cluster on a time-scale shorter than, or comparable to, the crossing 
time of stars in the cluster. This causes the cluster to expand, lose rnass or 



even completely dissol ve (e.g. iHilld Il980l : iGoodwiril Il997l : iKroupa et al.l 12001 



iGever Sz Burker"tll200ll ). The simulations of clusters including residual gas ex- 
pulsion generally assume a spherically symmetric gas configuration and instanta- 
neous gas removal, resulting in an unbound cluster whenever the star formation 
efficiency is less than 50%, i.e. when the total gas mass that is expelled equals 
the mass of the stars. Detailed hydrodynamic calculations, taking into account 
the dumpiness of the gas, show that only a small fraction of the gas is acceler- 
ated to high velocities, escaping through irregular outflow channels, causing a 
larger fraction of stars to remain bound ( Dale et al. 20051 and Clarke in these 
proceedings). From observations it appears that some clusters remove their pri- 
mordial gas more efficiently than others. For example, the extremely young 
(1 ± IMyr, IStolte et alll2004l ) Galactic cluster NGC 3603 has cleared its inner 
region, while NGC 346 in the Small Magella nic Cloud is still s urrounded by gas, 
even in its core, whereas its age is 3=b 1 Myr (jSabbi et al.ll2007l and Nota in these 
proceedings). 



Summary 
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Rapidly expanding clusters have a wide range of implications on the evo- 
lution of galactic discs (|KrouDal [20051 ^ . the CIMF (jBoilv & KrouDa"2003") and 
the early disruption, or infant mortality, of clusters (La da &: Lada 2003). In 
addition, estimates of dynamical masses at young ages (< 30 Myr) are affected, 
since these rely on the assumption that the cluster is in vir ial equilibrium, which 
is not the case right after gas expulsion (Goodwin &: Bastian 2006,). 

Several observational confirmations of the infant mortality scenario have 
become available recently. The most direct one comes fror n the age distr i bution 
(dN/dt) of (young) clusters. For the solar neig hbourhood fLada k Ladal (|2003l ) 
find that there is a steep drop in dN/dt going from the embedded (< 3Myr) 
clusters to the young clusters that have expelled their gas (> 3Myr), corre- 
sponding to ~ 90% infant mortality. This value was confirmed independently, 
by determining the mean formation rate of bound clusters in the last few Gyrs, 
correcting for various disruption processes, and comparing this to the total star 
formation rate following from field stars ( Lamers Sz Giele"i 2006 and in these 
proceedings). 

From the age distribu tion of (optically detected) clusters. iFall et all ( 20051 ) 
and Bastian et al. ( 2005al ) derive similar values for the infant mortality rates for 
the clusters in t he Antennae galaxie s and M51, respectively. A novel approach 
is presented bv iPellerin et al. ( 20071 ) (also in these proceedings). They study 
groupings of individual stars of different spectral type in nearby (few Mpc) spiral 
galaxies using HST/AGS data. The O stars are strongly clustered, while the B 
stars are more equally spread out over the disc, suggesting that infant mortality 
works on time-scales comparable to the life-tir ne of O stars. From the surface 
brightness profiles of slightly resolved clusters, Bastian fc GoodwinI (j2006) find 
evidence of escaping stars in the outer halo of young (few Myrs) clusters from 
high resolution HST/HRC imaging. 



4. The evolution of the (globular) cluster mass function 



Larsen (in these proceedings) highlights some of the standing problems in our un- 
derstanding of the globular cluster mass function (GCMF). The general picture 
is that clusters form from a universal power-law CIMF wi th index — o: = —2 , 
as is observed for young clusters in several galaxies (e.g. I Zhang k Falll 119991 : 



I _ ,1 1 - ^ 

de Griis et al.M2003l : iBik et al.M2003l ). Through dynamical evolution, the low- 



mass clusters are preferentially destroyed, causi ng the GCMF to turn-over (e.g. 
Chernoff fc Weinbergl Il990l : iFall fc Zhang] l200ll ) . Although the general idea is 
quite established, there are still several issues to be resolved. First, if dynamical 
evolution removes low-mass clusters, then the location of the turn-over in the 
GCMF, i.e. the turn-over mass, will depend on the strength of the tidal field. 
This implies a decrease of the turn-over mass with galactocentric distance (Rq), 
because the strength of the tidal field decreases with Rq. 

However, from observations it follows that the turn-over mass i s almost 
unive rsal among different galaxy types and independent of Rq (e.g. iRichtleil 
20031 ). Solutions to overcome this problem have been introduced in analyti- 
cal models, such as including a strongly increasing velocity anisotropy with Rq 
( Fall k Zhang 200ll ). Observations of globular cluster kinematics have shown 
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that these assumptions are not reahstic (e.g. ICote et al.ll200ll for the M87 clus- 
ters). 

The term "turn-over" can be misleading. The dA^/dlogM distribution, or 
the magnitude distribution (since the magnitude scales with the logarithm of the 
mass, assuming a constant M/L-ratio), is peaked. A dN/dM representation of 
the GCMF is flat on the low-mass end and on the high-mass end it approaches 
the (initial?) power-law distribution with index —2 (see Fig. 2 of Larsen in these 
proceedings). 

If the disruption time of clusters (tdis) depends on their mass as a power- 
law with index 7, then the low-mass end of a sin gle-age population evolves to a 
power-law with index 1 — 7 ( Lamers et al.|[2b05al ). A flat dN/dM distribution is 
thus achieved when 7 = 1. It is worth noting that the final index is independent 
of the initial index a, implying that the low-mass end of the globular cluster MF 
reflects only disruption and not formation. This holds of course only if the initial 
index —a is smaller than 1 — 7. T h e evo lution of GCMFs that are bell-shaped 
initially are discussed bv lVesperinil ( 20001 ). 

For a multi-age population which has formed with a continuous rate and 
for which the age spread is much larger than the typical tnisi the index of the 
dN/dM distribution approaches a value of 7 — a (iBoutloukos fc Lamersl[2003l ) 
(BL03), which is —1 in the case of —a = —2 and 7=1, i.e. steeper by an index 
1 as compared to the single-age case. This type of multi-age populations were 
studied by BL03 and they find a mean index of the mass function of ~ —1.4, 
implying 7 ~ 0.6, for the cluster populations in the SMC, M33, M51 and the 
solar neighbourhood. 



5. Some notes on 7 



The value of 7 = 1 is needed to get the flat dN/dM distribution of the GCMF 
on the low- mass end, while BL03 find 7 ~ 0.6 from the dN/dM distribution of 
clusters with ages up to a few Gyrs. Did globular clusters evolve differently than 
their younger counterparts? Or did globular clusters not form with the same 
initial mass function as clusters now? 

Th ere is theore t ical s upport for a value of 7 = 0.6. The arguments are as 
follows: iBaumgardtl (|200ll ) showed that tjis depends on the relaxation time (t^h) 
and the crossing time (tcr) of stars in the cluster as tdis — ^c^^ which, with 



the expressions for t^h and tcr from ISpitzerl ( 19871 ). can be reduced to 



tdis OC 



N ' 



M 



1/2 



(1) 



with A the Coulomb logarithm which scales with as A :^ O.IN ( Spitzei^' l987l ) 
and Th the half-mass radius of the cluster. iBaumgardt &: Making (2003 ) veri- 
fied this relation using detailed A-body simulations of clusters dissolving in the 
Galactic tidal field and they found that x ~ 0.75. They assumed that clusters 
are initially in tidal equilibrium with the Galaxy, implying a constant cluster 
density, i.e. rh oc M'^ with A = 1/3, at a given galactocentric distance. This 
re duces Eq.fflto tnis oc [A/In A ] The sam e scalin g of t^is with N was found 



bv IVesperini &: Heggid ( 19971 ). iLamers et al.l ( 2005bl ) showed that this relation 
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can be well approximated by tjis oc A^'''^, agreeing with the empirical findings 
of BL03. The scaling of A with N is often ignored by taking In A constant (e.g. 
Fall &: Zhand i2001). Combined with the assumption of a constant cluster den- 



sity, i.e. A = 1/3, then this results in a linear scaling of tdis with N (Eq. [T]) 
equivalen t to 7 = 1 , which is required to model the low-mass end of the GCMF 
(jWaters et al .1120061 ). However, this value of 7 does not follow from recent A^- 
body simulations, but from the assumptions that t^is scales linearly with ifh, a 

contant A and a consta nt cluster density. 

The assumption of iBaumgardt &: Makind ( 20031 ) and IVesperini Sz Heggie 



( I997I ) that clusters start tidally limited, i.e. A = 1/3, is also not realistic 



From recent observations it follows that the relation between rh and Af is muc h 
shallower, i.e. A ~ 0.1 (e.g. IZeof et al.l[T999l : lLarsenll2004l : iBastian et al.l[2005bl l. 



With such a mass-radius relation, the massive clusters ar e not filling their "Roche 
lobe" . This situation was considered in A^-simulations bv lTanikawa Sz Fukushig^ 
( 2OO5I ) . They modeled the disruption of clusters in weak tidal fields, where the 



tidal radius of the clusters is smaller than the Roche lobe radius. They find that 
in those 0.9 — 1. 0, i.e. almost as if the cluster is evolving in isolation 

(then X = 1, ISDitzed[l987l l. Hence, x depen ds on A . Lets consider the case of 



a realistic mass-radius relation, i.e. A = 0.1 ( Larsen| [2004). Assume a constant 



mean stellar mass, i.e. M ex N and lets approximate A^/lnA by N^-^^ ^ which is 
a good approximation for N » N/A ~ 10. This, combined with x ~ 1 and 
Eq. [T] results in 7 ~ 0.5. When clusters are disrupted by external perturbations 
(e.g. disc/bulge shocks, encounters with giant molecular clouds or spiral arms) , 
tdis scales with th e cluster density (pc) (e.g. ISpitzeij[l987l : [Cieles et aLll2006d : 



iGieles et al.ll2007l ). which combined with A = 0.1 results in tdis oc Ad ' , again a 
value for 7 that is smaller than 1. 

In neither of the cases described above we get 7 = 1 needed for the low-mass 
end of the GCMF. Hence not only the universality of the turn-over mass, but 
also the shape of the GCMF can not be explained by the standard picture where 
clusters start with a power-law CIMF with index —2 and evolve dynamically. 
The three predicted values for 7 are quite close to the observationally determined 
value of 7 = 0.62 ± 0.06 (BL03). 
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